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A trajectory of a harmonic oscillator obeying the Schrodinger wave equation is 
exactly derived and illustrated. The trajectory resembles well the classical orbit be- 
tween the turning points, and also runs through the tunneling region. The dynamics 
of the 'particle' motion and the wave function associated with the motion are pro- 
posed. The period of a round trip on the trajectory is exactly equal to that obtained 
in classical mechanics. 



1 Introduction 

In general the motion of a particle whose beam shows an interference phenomenon is de- 
scribed by the wave equation. On the other hand the orbit of the particle described in 
classical mechanics plays an important role in the fields such as electron optics, particle 
accelerators, radiation from electron beams and so forth. The significance of the scattering 
wave function, wave-optical approach, was clarified by comparison with the classical wave 
function, ray-optical approach, by Gordon. [|I|] 

Geometrical optics describes the movement of the light corpuscle in the space as the ray 
of light. Formulas for the ray of light derive from the eikonal equation. This corresponds 
to the Hamilton characteristic function for the particle motion, which can be derived 
from the Schrodinger wave function in the WKB approximation for the system with the 
stationary potential. 

The concept of the orbit or the ray of light in the wave phenomena gives a comprehensive 
image and physical insight of the process, although it has been an approximate idea. ^, |^ 
The field ion microscope and scanning tunneling microscope are now available to image an 
atomic movement on solid surfaces. @, An atom can be picked up and moved to an 
arbitrary place. These suggest the possibility of the precise description of the motion of 
the particle even in the atomic scale. 
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A proposal of extension of the light ray to the shadow region had been made in order to 
clarify the physics of the diffraction. The relation between the wave and the light ray or 
orbit should be investigated more carefully. Trajectory in the optical wave, extension of a 
ray of light, was developed by generalizing the eikonal to the mode characteristic function 
with dynamical assumptions. ||n| 

If the motion of the particle is restricted to that in the classical region where classical 
mechanics is the case, the orbital motion is derived from the Hamilton characteristic func- 



tion by the Hamilton- Jacobi theory for the system with the stationary potential. |12] If the 
characteristic function could be generalized to that valid for every space region, it would 
be able to pinpoint the motion of the particle even in the tunneling region. [|l^ 

Along the similar way of thinking, Bohm proposed a quantum theory with "hidden 
variables" to suggest objective description of individual systems at a quantum level of 
accuracy. [14| It has made a general scheme of the causal interpretation but not given a 



concrete trajectory with the dynamical behavior of a 'particle' in the space-time region from 
the wave equation. The goal of every causal theory or interpretation as to the quantum 
mechanics is a complete description of an individual real situation as it exists independently 
of acts of observation. |jl5| The motivation is summarized in the Einstein's feeling that the 
statistical prediction of the quantum theory is correct but by supplying the missing elements, 
it could be in principle got beyond statistics to a determinate theory. [0 To fulfil the aim 



it might be necessary to exploit the suitable mathematical tools to describe the particle 
motion in quantum mechanics. 

In the present paper, a trajectory of the harmonic oscillator in one and two dimension 
is exactly derived from the Schrodinger wave function. The trajectory and dynamical 
motion are compared with the classical orbit and dynamical behavior in the classical region. 
The relation between the traveling waves associated with the 'particle' motion and usual 
stationary wave function is discussed. 



2 Dynamics and wave function 

The way to derive the trajectory from the wave equation is described and the significance of 
the wave function is discussed. The Schrodinger wave equation for a particle in a stationary 
potential, V, in one spatial dimension, is 

The equation is separable in t and x. A wave function for a 'stationary' state is written as 

= e^^{-iEt/h)^{x,E), (2) 

where is a constant of separation. Function $ is a general solution of an ordinary differ- 
ential equation of second order. The boundary condition that the function be bounded and 
continuous everywhere in the domain defined is not imposed on the function yet. Function 
$ is the 'eigenfunction' of the Hamiltonian operator H . Constant E is the 'eigenvalue' and 
the energy of the 'stationary' state. The wave function ^> determines the 'eigenstate', or 
the 'mode', specified by E. The constant E should be called the mode parameter. 
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If the parameter takes specific, discrete, values, the wave function satisfies the boundary 
condition. The state with a specific parameter is usually called a stationary state or the 
mode, and function $ is the eigenfunction. For the present, the parameter is assumed to 
be a real number and the words, such as stationary, eigen, mode and particle are used with 
quotation marks. 

The wave function of the form of 



E) = ^p{x, E) exp[i W{x, E)] (3) 
is sought, where functions p and W are real number. Function W satisfies approximately 

in the WKB approximation. The region where this equation holds should be called the 



classical region. This is the Hamilton- Jacobi equation with energy E. [|T^] If the solution 
for hW of this equation is written as Wd, the Hamilton characteristic function is given 

by m 

Wci{x,E) = ^2m{E - V)dx. (5) 



Function W satisfies the nonlinear equation (4) in the reference of Bohm |T^, expressed in 



the form S = —Et+hW, the Hamilton- Jacobi-like equation with the "quantum-mechanical" 
potential. Function W is assumed to satisfy the condition that in the classical region of x 

hW{x,E) c:^W{x,E)d, (6) 

except an additional constant independent of x and the mode parameter. If function W 
can be determined uniquely, function W should be named the mode characteristic function 
(mcf) for the system. The mcf W{x) is derived from the wave function as follows. 
Since the wave equation is the ordinary differential equation of 2nd order, the function $ 
is composed of two linearly independent solutions, say ui and U2, 

^{x,E) = aui{x,E) +bu2{x,E), (7) 

where a and b are complex constants. Constants a and b should be determined by the 
following assumptions from the theoretical point of view. 

Assumption 1: Function hW in the classical region must be as approximate to the 
corresponding characteristic function Wd as possible in the sense of relation (^. 

Assumption 2: The results derived from the equations of motion for W defined by Eq. 

should be as approximate to the ones in classical mechanics as possible. 

Function W{x) should be determined as the phase or argument of function $ like Eq.(^. 
The function \l/ consisting of this $ may be said to represent the traveling wave associated 
with the motion of a 'particle' in the 'mode' E. 

The equation of motion for the 'particle' is assumed 

dW 1, , 
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where to is a constant (independent of t) that should be determined by the initial condition 
for the system. Variable t should be the dynamical time for the system. 

Dynamical time t should be determined so as to increase monotonically as the 'particle' 
moves. From equation (|]), a phase velocity 

dt [ dxdE ) ' 

is obtained. If the 'particle' starts from a position, say xq, at an initial time and dt/dx > 
for X < Xb, it could be considered that the 'particle' runs from xq to Xb- Here Xb stands 
for a turning point or an endpoint of the potential. After reaching Xb, the 'particle' returns 
to Xq with the mcf of —W{x, E) + 2W{xb, E), which guarantees the monotone increase of 
time t given by equation (||). 

The wave function associated with the 'particle' motion in a bound state should be 
described as follows. If the function $ (|^) represents the motion of the 'particle' traveling 
to the right, function $* = y/pexp[—iW], which is also the solution of the wave equation, 
stands for the motion to the reverse direction. Let the endpoints of the potential in the x 
coordinate be a and b. Then the 'particle' moves in the region a < x < b. Let it start from 
a point Xq to the positive x direction and the mcf be W{x). The traveling wave associated 



with the returning motion from b to Xq or a should be given by y p(x, E) exp{i [—W{x, E) + 
2W{b,E)]). 

The wave function observed at x(> xq) is assumed to be the superposition of the trav- 
eling waves of either motion 

^{x,E) -^*{x,E)exp[i2W{b,E)]. (10) 

This is finite, zero, at b. If the 'particle' turns at a and runs to Xq or b, the mcf should 
be given by W{x, E) — 2W{a, E) + 2W{b, E). The wave function for x(< Xq) is thus to be 
written as 

- $*(a;, E) exp[i 2W{b, E)] + $(x, E) exp {i [-2W{a, E) + 2W{b, E)]) . (11) 
This is finite, zero, at a. The wave functions (|10|) and (111]) are bounded for a < x < 6 for 



any 'mode', although they might not be continuous at Xq. 

For X (> Xq) the traveling wave associated with a round trip of the 'particle' gets the 
shift in phase by —2W{a, E) + 2W{b, E). The wave function associated with round trips 
of the 'particle' motion would result in, being averaged for one cycle, 

1 l-exp{tN[-2W{a,E) + 2W{b,E)]) 
N l-exp{i[-2W{a,E) +2W{b,E)]) 

X {^{x,E)-^*{x,E)exp[i2W{b,E)]). (12) 

If the number N is large, this shows a sharp resonance if E satisfies 

- 2W{a, E) + 2W{b, E) = 27c x integer. (13) 

This resonance condition gives rise to the stationary state and the eigenfunction. This is 
the exact version of the Bohr-Sommerfeld quantum condition. An approximate but a little 
general expression for it had been presented previously. |]17| It might as well be interpreted 



that the observed wave should be proportional to the wave function mentioned above. 
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3 Simple harmonic oscillator 



The trajectory of the simple harmonic oscillator is discussed. The wave equation for the 
oscillator with mass m for energy E is written as 



+ 1 ^(x) = E^(x) 



\ 2m dx"^ 2 

A solution of this equation of the form (|^ that leads to function (|^) is 



(14) 



$(+), 



x] 



M(a, 1/2, z) - I , , V(a, V2, x) 



v^r(V2 - a 
where T{x) is the gamma function and 

a = 1/4 (1 — 2E/huj) , z = mujx'^ /fi 



exp(-2;/2), 



(15) 



(16) 



Functions M and V are linearly independent confluent hypergeometric functions. |18| The 
former is the Kummer function, which is expressed in the series form 



M(a,V2,^) = E^ 



(V2)fcA:! 



;i7) 



and the latter is defined as fill 



V^(a,V2,x) = T{-^/2)^muj/hxM{^/2 + 0,3/2,^). (18) 

Expression ([T5|) is the 'eigenf unction' with 'eigenvalue' E, not always bounded at x = cxd, 
of Eq. (|ll. 

The mcf W{x, E) satisfying two assumptions mentioned in section is given by 



W{x, E) = arctan 



-r(l-a) V{a,%x) 
v/^r(V2-a) M{a,% z) 
arctan[F(a, 1/2, x)] 



(19) 



This mcf multiplied by % is well approximate to Wc\{x, E) with iyci(0, -E) = in the classical 
region or in the neighborhood of x = 0. The validity of the mcf will be recognized in the 
following discussion. 

The equation of motion is given by Eq. (^, or 



1 

X 



■v^(l - a) - nh - a) - ^^^g^(V^ + "'^/^'^) + ^l°g^(-'V2, 



da 



da 



rfl - a) 



M(V2 + 0,3/2,^) 



r(V2 - a) M(a,V2,^)[l + F(a,V2,x) 



(20) 



where \E'(x) = d\ogV{x) / dx is the psi function. [|T^ The oscillator has been assumed to be 
at X = at t = 0. 
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By using the asymptotic form of the confluent hypergeometric functions, it is ob- 
tained for X large 



F(a,y2, x) ^ tan 



TT / 2E\ 



(21) 



It is thus obtained that t{x = oo) = tt/2uj. The x dependence of function t{x) is monotone 
everywhere, as proved by a computer calculation. Therefore it can be considered that the 
oscillator moves between the end points, x = — oo and oo, without interruption. 

If the oscillator starts from a; = at t = and goes to x = oo, it returns there and 
comes back to x = — oo. Then it goes back to x = 0. The mcf W{x, E) for one cycle is 
written as follows: 

{arctan[F(a,y2, x)] (x = ^ oo), 

2iy(oo, E) - arctan[F(a, 1/2, x)] (x = 00 ^ -00), (22) 
AW{oo,E) + arctan[F(a,V2,x)] (x = -00 ^ x = 0). 

These are determined so that time t increases as the oscillator runs along the trajectory. 
From the asymptotic form of ^(0,^/2, x) for x large, (PH), and ^(0,1/2, 0) = 0, it is seen that 
the period of the motion is equal to 27r/uj that is just equal to that in classical mechanics. 
The oscillator runs throughout the space between x = — 00 and 00 and the phase velocity 
dx/dt in the tunneling region is very large. 

The variation of the mcf after one cycle stands for the change of the phase of the wave 
function. It is at any point 

- 2W{-oo, E) + 2W{oo, E) = 4W{oo, E) = n {1 + 2E/nuj) . (23) 

The resonance condition of the wave function (0) leads to the eigenvalues of energy. Thus 
it holds for the stationary state 

E = /icj (n + 1/2) , n = 0, 1,2,---. (24) 

If equation ( PU] ) is solved inversely for x as a function of t, the function x(t) is well 
approximated by the classical oscillation 



X = ^2E/mu^ smut. (25) 

Expressions ( pOD and (|25D are illustrated in Fig. |l|. The discrepancy between the two occurs 
at about times, t = 7t/2ujx half integer, when the oscillator is running through the tunneling 
region. 

The similarity suggests that the mcf (0), satisfying assumptions in section 0, is the 
correct one for the harmonic oscillator. 



4 Wave function 

The wave function (|15[) multiplied by exp(— lEt/Ti) describes a wave traveling to the right. 
According to the mcf (|22D it may be considered that it represents the wave associated with 
the oscillator motion for x = to 00. 
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Figure 1: Position versus time of a harmonic oscillator with E = ^/2fiuj. Solid line stands 
for the trajectory and broken line for the corresponding classical orbit. 

A wave traveling to the left associated with the oscillator motion from x = cxd to — oo 
should be given by 



M{a,y2,z) + i 



ra-a) 



-r(a,V2,x) 



v^r(V2 - a) 
X exp[-z/2 + i2W{oo,E)]. 

A linear combination with ([T5|) like expression (|IOD, 

= 2 exp(— ma — 1/22;) M(a, 1/2, z) cos vra 

rfi - a) 



(26) 



+ 



v^r(V2-a 



■V"(a,y2, x) sinvra 



(27) 



gives rise to a wave function finite for < x < 00. The boundedness at a; = 00 is seen as 
follows. Expression ( pT]) can be rewritten as 



^i{x,E) 



20F 



r(V2 - a) 



exp(— i na — z/2)xU (a + ^2,^2, z). 



{2i 



large []TB[ it is found that 



Function U is the Kummer function. |18] From the asymptotic form of function U for x 

uhat 

$i(x, E) ~ ^^f" exp(-i vra - 1/2^), (29) 



r(V2-a)- 

which tends to zero as x tends to infinity. 

For X negative, an associated wave for the oscillator motion reflected at x = —00 and 
going back to x = should be written as 



$(+)(x,E)' 



M(a,V2,^) 



r(i - a) 
Ar(V2 - a, 



■V(a,V2,x) 



exp[-z/2 + i4Vr(oo,^)]. (30) 



The wave function for x negative is 

$2(a;, E) = {x, E)' - {x, E) 
= 2 exp(— zSvra — 1/22;) M(a, 1/2, -z) cos vra 
Til -a) 



+ 



0Fr(V2-a 



-V"(a,y2, x) sin vra 



(31) 



Since M(a,y2,2;) is an even function of x and ^(0,1/2, x) is an odd, there is a symmetry 
between functions $1 and $2 



^2{x,E) = ^i{-x,E) exp{-i2na). 



(32) 



A wave function $1 except a constant but ii^-dependent factor is shown in Fig. |^ for several 
E's. 




Figure 2: Wave function $1(0;, i?)e*™ vs x for E = i^/icjx (a) 1, (b)1.5, (c) 2, (d) 2.5, 
(e) 3. 

It can be seen from the figure that function $1 is smoothly continuous to function $2 
at a: = if £^ = 1/2/10; x [1 or 3], which is the eigenstate. 

If 4W{oo, E) is a multiple of 27r, in general, expression ( ^Tf ) becomes smoothly contin- 
uous to equation (^) at a; = 0. The two expressions constitute the eigenfunction for the 
stationary state finite and smoothly continuous for any space point with the eigenvalue E 
discrete. 



5 Two dimensional motion 

The simple harmonic oscillator runs only between the end points, x = —00 and x = 00. 
Except that it runs through tunneling regions, it oscillates like a classical particle. Here, 
the harmonic oscillator moving in the space of two dimension is studied. The extension is 
straightforward if the partial differential equation is separable in variables. 
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The wave equation for the harmonic oscillator in two dimensional space with energy E 
is given by 



2m ydx"^ dy 



Hx,y), (33) 



where m is the mass, and ooi and UJ2 are proper frequencies. By introducing a constant of 
separation of variables, E2, the above equation can be decomposed into two equations 

Ei$i(x,i?i) = f-_^ + ^a;2^M^i(x,i?0, (34) 

E2^2[y.E2) = (-^^ + Y^2Vj$2(2/,i?2), (35) 

where Ei = E — E2. 

By following the discussion of section ^, the mcf in the 'mode' {E, E2), is obtained as 

W{x, y, E, E2) = Wi{x, E, E2) + W2{y, ^2), (36) 

where 

Wi{x,E,E2) = arctanF (1/4(1 - 2Ei/^cui), 1/2, x) , (37) 

W2{y, E2) = arctanF (1/4(1 - 2E2/huj2),y2, y) ■ (38) 

The equations of motion are given by equation (H) , or 

,^,^^, mm) ^,^^,SW^^ (39) 

where to and t2 are constants to be determined by the initial condition. The whole trajectory 
for —00 < X < 00 and —00 < y < 00 can be obtained by using mcf's (|22| ) for x and y 
coordinates. The projected motion of the 'particle' onto the x or y coordinate is the 
periodic one with period of 2tt/ui or 2tt/uj2, respectively. If the ratio 001/002 is irrational, 
the trajectory in the two dimensional space is not closed as is the case in classical mechanics. 

An example of the trajectory for the system with parameters Ei = 3/2^1^1 and E2 = 2huj2 
and 001/002 = 0.8 is shown in Fig. ^. The oscillator is set to start from the origin at t = 0, 
or to = ^2 = 0. In the figure the corresponding classical orbit with the same parameters is 
also drawn. It could be seen that the bigger the values of the parameters, the closer the 
trajectory and the classical orbit with each other, which shows the correspondence principle. 



6 Conclusion 



By generalizing the argument on optical wave, |Tl]| the dynamics that should figure out a ray 



of 'particle', a trajectory, in a 'mode' of the Schrodinger wave equation of the completely 
separable form, especially in one spatial dimension, has been proposed. The dynamics 
should be determined by the mcf and dynamical assumptions on it. The trajectory thus 
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Figure 3: The trajectory of a two-dimensional harmonic oscillator with Ei = ^/2huJi, 
E2 = V2^i^2 and uoxjuJi = 0.8 (solid line), and the corresponding classical orbit (broken 
line) . 

determined resembles well the orbit of the corresponding state in classical mechanics in the 
classical region. It runs also through a tunneling region. This is verified for the harmonic 
oscillator system. The period of a 'particle' in the oscillator is exactly equal to that in 
classical mechanics. 

The wave function bounded everywhere but not always continuous for any bound state 
could be made by superposing the traveling waves associated with the 'particle' motion. 
The eigenfunction of the stationary state should be interpreted to be the wave function of 
the resonating state in the potential of the system. It is to be noted that all solutions of the 
wave equation are necessary in order to derive the mcf to get the trajectory. This suggests 
the role of all the solutions of the wave equation. Since the trajectory is determined by 
the phase or argument of the traveling wave function, it does not always show all the path 
of the energy transfer, which might be necessary for the real particle motion in the wave 
equation. 

The accordance of these characteristics between classical and quantum mechanics sug- 
gests the validity of the dynamics defined here. The main difference from the Bohm's 
theory [|I^ stems from the dynamical assumption (^. The significance and the consistency 



with the quantum theory should be verified for a more important system such as the one 
with the Coulomb potential. By treating the scattering problem, the statistical but an in 
principle determinate nature in quantum theory will be shown on getting the cross section. 

It suggests the consistent existence of the trajectory in wave mechanics and the 
significance of the relation between the particle motion and the wave function, although 
there may remain a lot to be considered about the observation process of particle and wave 
phenomena. 
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